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Optical imaging plays a critical role in advancing our understanding of three dimensional dy¬ 
namics of biological systems. Coherent imaging (Cl) methods exploit spatial phase information, 
encoded through propagation of coherent signal light emerging from a specimen, to extract a three- 
dimensional representation of the object from a single high-speed measurement. Until now, Cl 
methods could not be applied to incoherent light, severely limiting their ability to image the most 
powerful biological probes available - fluorescent molecules - with sufficient speed and volume to ob¬ 
serve important processes, such as neural processing in live specimens. We introduce a new imaging 
technique that transfers the spatial propagation phase of coherent illumination light to incoherent 
fluorescent light emission. The transfer of propagation phase allows Cl techniques to be applied to 
fluorescent light imaging, and leads to large increases in imaging speed and depth of Held. With 
this advance, biological imaging of fluorescent molecules is significantly expanded. 


INTRODUCTION 

The ability to spatially resolve the structure of objects plays a critical role in understanding nearly every physical 
system, providing insights into fields as diverse as astrophysics mm , molecular solvation dynamics [3j, geophysical 
systems |f], the behavior of solid-state materials [5], and biological systems [5]. Likewise, the ability to resolve 
temporal dynamics is essential for unravelling processes of physical systems in relation to internal and external stimuli. 
Obtaining a comprehensive understanding of the interplay between spatial organization and functional behavior in 
dynamic specimens remains a pressing challenge, particularly in biological systems. 

Many biological processes are mediated chemically, and relevant mechanisms are often revealed only by labelling 
such processes with fluorescent probes [7]. Biological sciences rely on the molecular specificity of fluorescent probes 
to study targeted chemical processes. Three dimensional images formed by collection of incoherent light can be 
obtained with high spatial resolution using techniques that restrict illumination and/or detection to a small volume 
in the object region [HI IS]. However, restricting data collection to a small volume necessitates serial data acquisition, 
limiting imaging speed and degrading the ability to track dynamic behaviors in 3D. 

With notable exception [10], coherent radiation enables much faster 3D imaging than is currently possible with 
incoherent light. Three dimensional images of an object can be reconstructed from a single, rapidly acquired, two 
dimensional image of light collected from spatial points throughout the object by exploiting our understanding of 
coherent wave diffraction behavior. This powerful tool has proven useful in seismic, ultrasonic m , photo-acoustic 
m , optical nanu, x-ray Hang, and electron-beam m imaging systems for rapid detection of spatial organization 
and temporal dynamics. These techniques, generally referred to as coherent imaging (Cl), rely on the recording, or 
numerical recovery, of optical spatial phase information that reveals the distance the light has travelled from the point 
at which it scatters in an object to where it is recorded. 

Special techniques are required to extract the phase and amplitude of the diffracted signal field because optical 
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detectors are incapable of responding to the phase of a light beam. Optical interference converts phase differences 
between two beams to intensity variations that can be recorded and used to recover the amplitude and phase. Holo¬ 
graphic microscopy (HM) pL3] represents the current state-of-the-art [18]. Alternatively, the amplitude and phase of 
a coherent field diffracted from an object can be recovered through inversion of coherent light propagation with a 
regularized, iterative optimization algorithm mm, without requiring interference. 

Despite the advances of these techniques for imaging coherent radiation in a three dimensional region, the problem 
of dynamic volumetric imaging of 3D incoherent light emission, such as from fluorophores, remains vexing. The phase 
of the fluorescent light emitted by the molecules bears no stable relationship to the illumination light or neighboring 
fluorophores in the specimen, rendering the fluorescent light completely spatially incoherent. As a result, no phase 
information is available to pinpoint the location of the fluorescent emitter, normally excluding the application of Cl 
techniques to fluorescent imaging. The vast potential of Cl has lead to its application to holography using incoherent 
sources [lOl - E!?| and fluorescent emission j23fEH] . Such techniques fail to image in the presence of optical scattering, 
with fluorescent light emission at high speed, or with large volumes in an epi-collection configuration, all of which are 
needed for monitoring dynamic biological behaviors. 

Cl is able to image 3D volumes at high speeds, but these powerful tools have not been applicable to incoherent 
light imaging until now. In this Article, we introduce an imaging method that applies the high speed 3D imaging 
capabilities of Cl to incoherent fluorescent light emission. Encoding the spatial phase of a coherent illumination source 
into a measurement of fluorescent light intensity collected with a single-element photodetector makes this possible. 
Dubbed Coherent Holographic Image Reconstruction by Phase Transfer (CHIRPT), this technique enables application 
of the full range of strategies used for coherent volume imaging with incoherent light. CHIRPT can, in principle, be 
applied to any incoherent emission or scattering process that is illuminated by a coherent source. CHIRPT imaging 
will prove particularly valuable for biological applications because of its ability to image at high acquisition rates in 
an epi-fluorescence configuration, its relative insensitivity to scattering due to single-element detection, and ability to 
image a large depth-of-field. 


PHYSICAL PRINCIPLES OF CHIRPT 

CHIRPT enables Cl techniques to be used with incoherent light. The illuminating light is modulated with a temporal 
intensity pattern unique to each spatial point in the object region. The modulation is formed through interference 
of coherent illumination beams, and this interference transfers the spatial propagation phase of the illuminating 
fields to the temporal modulation pattern. The modulated illumination light excites fluorescent molecules, and the 
light emitted by those molecules mirrors the temporal illumination intensity pattern at the molecule location. The 
fluorescent light emitted by the object is collected onto a single element photodetector, and the signal generated by 
the photodetector contains information that identifies the position of fluorescent molecules in the object. The signal 
is a compilation of intensity modulation patterns emitted by fluorescent molecules in the object. With this set of 
modulation patterns, spatial propagation phase in the object can be recovered through a simple, direct computation. 
Using this spatial phase, we can now exploit the power of Cl techniques to form images with incoherent light. 

Generating the intensity modulation pattern is the key aspect of CHIRPT, which is shown conceptually in Fig. [T] 
We begin by directing a laser beam into a cylindrical lens to form a line focus. The line focus passes through a 
modulator disk spun at a constant angular velocity. The modulator disk is lithographically printed with a mask that 
contains a transmission grating designed to vary the groove density (spatial frequency) with angular position (Fig. SI). 
As the laser line passes through the spinning disk, a portion of the laser power diffracts to an angle, determined 

by the groove density. The undiffracted zero order beam is undeviated upon propagation through the mask, and thus 
stationary through the duration of the scan. The two beams - diffracted and undiffracted - are imaged to the object 
region as light sheets whose crossing angle varies with time, as represented in Fig. UK It is the interference between 
these beams that creates the intensity modulation pattern. 

The evolution of the interference pattern in the object plane ( x , z) with time is dictated by the relative optical 
path, AOPL(x, z; t), accumulated by the two beams as they propagate from the mask plane to the object region. We 
can understand the OPL acquired by each beam using the framework provided by Fermat’s principle of least action 
)2~6] . The OPL of the undiffracted beam is fixed, whereas the diffracted beam OPL depends on the diffraction angle. 
An ideal imaging system introduces zero OPL difference for paths taken between conjugate points in the modulation 
and object planes. OPL variations caused by imperfect imaging are captured by the pupil aberration phase, </> pU pii(0- 
This aberration phase causes further variation of the illumination intensity, which can be numerically corrected with 
CHIRPT (Fig. S4). 

Along with the imaging system OPL, the beams pick up additional OPL as they diffract along the optic axis ( z ) away 
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FIG. 1: Principle of CHIRPT imaging, (a) The angle of diffraction from the modulator, d\(t), varies with scan 
time with periodicity T = \/v r . Scan times of —T/4 and Tj2 are indicated by the red and blue beams, respectively. 

(b) The difference in optical path length, AOPL/A, between the diffracted and undiffracted beams at t = —T/4, and 

(c) t = T/2. Contour lines represent integer values of AOPL/A, corresponding to peaks in the illumination intensity 
pattern, (d) AOPL/A as a function of scan time for the points (pi,P 2 ,P 3 ) plotted in (b) and (c). Here, AOPL was 

computed assuming a paraxial imaging system, where A4>(x, z ; t) « 2n Ax (M Kit) — it A A z (M k\ t) 2 . 


from the object image plane. The contribution of the imaging system OPL and the diffraction OPL determines the total 
OPL difference, AOPL, between the two beams. A derivation of the illumination intensity modulations, Im(x,z;t), 
caused by AOPL shows the intensity follows the form Im(x,Z',t) oc cos[A4>(a;, z; f)]. Here, the phase difference is a 
function of AOPL and the wavelength of the illumination light, A, given by A 4 ){x,z\t) = 2tt AOPL(:r, z\ t)/A. 

With an understanding of the contributions to the OPL described above, we can now write out a simple model of 
the CHIRPT imaging process. The relative phase is given by (Supplementary Information): 


A4>(:r, z\ t ) = A 4> a (t) + 2ir Ax ( —M Kit) + 2-7T —- A z 

J A 

where A4> a (f) is the difference in aberration phase between the undiffracted and scanned beams as a function of 
time, and n 2 and ni are the refractive indices in the object region and the mask region respectively. Here Ax is the 
lateral location relative to the centroid of the illumination intensity distribution; Ki = Afc v r is the chirp parameter 
introduced by the modulation mask, describing the linear change in temporal modulation frequency as a function 
of lateral coordinate; A k is a parameter of the mask, setting the highest spatial frequency of the pattern; v r is the 
rotational frequency of the disk; M is the magnification of the image relay optics, defined as the ratio tube lens focal 
length to the objective lens focal length; and A z is the axial location relative to the focal plane of the imaging system. 

The illumination intensity is composed of a temporally-varying component, I i(x, z; t), and a stationary component, 
/o (Supplementary Information). Only the time dependent portion of the illumination intensity forms CHIRPT 


1 1 — ( — M Ki t ) — 1 

ni 


(1) 
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images. Retaining just the time dependent portion leads to an illumination intensity of: 


Ii(x, z;t) = — cos < A<!> a (f) + 27 t Ax ( —M m t ) + 2-7T ^ Az 
TT 1 Vri! J A 

Equation Q describes how each position ( x , z) in the object plane experiences a unique illumination intensity 
modulation pattern that varies with scan time. At each scan time, the intensity variations show uniform, tilted 
fringes representing a single spatial frequency in the illumination intensity. The uniform fringe spacing is caused by 
a uniform variation in AOPL. Figure [T] > and Fig. JTJ- demonstrate this for the two scan times indicated in Fig. 
where AOPL was computed from the theoretical analysis. Contour lines on the OPL maps in Fig. [TJ> and Fig. [TJ 
correspond to the peaks of the illumination intensity pattern in the object region. 

The OPL variation evolves with time as the disk spins. Each point in the object region experiences a unique 
temporal variation in AOPL. A visual representation of this can be seen for the three points indicated in Fig. & 
and Fig. [TJ:, where the full temporal trace gives a variation used to uniquely identify each spatial location. These 
OPL variations drive the intensity modulations at each point. Light emitted by fluorescent molecules excited at each 
point follows the illumination intensity temporal variation, provided the modulation frequencies are below the inverse 
fluorescent lifetime (Supplementary Information). 

Most of the energy diffracted from the mask is contained in the undeviated zero-order beam (j = 0) and the first 
diffracted orders (j = ±1). If the pupil phase is symmetric about the optic axis, e.g., spherical aberration, the two 
diffracted orders carry redundant and time-reversed spatial phase information. This redundancy in phase information 
introduces a spatial phase sign ambiguity akin to the twin-image problem that arises in on-axis holographic imaging. 
Unambiguous recovery of the spatial phase of the illumination system requires retaining only one diffracted order in 
CHIRPT imaging. We accomplish this by obstructing one of the first diffracted orders, in practice the j = — 1 beam, 
in a Fourier plane of the modulation mask (Fig. [2]). A portion of the j = 0 and j = +1 beams are spatially filtered in 
a plane conjugate to the back focal plane of the objective lens. The modulation mask facilitates spatial filtering of one 
diffracted order because it imparts a constant, non-zero spatial frequency in the vertical dimension, allowing complete 
spatial separation of the positive and negative diffracted orders in the filtering plane (Supplementary Information). 




FIG. 2: Schematic of the CHIRPT setup, (a) A constant diffraction angle in the vertical dimension results in 
vertical offset of the diffracted beam from the undiffracted beam (Supplementary Information). A horizontal slit is 
placed near the pupil plane of the objective lens to select a portion of the undiffracted beam and one diffracted 
order, (b) The illumination objective is oriented such that the illumination beams propagate along y = 0. (c) The 
diffracted and undiffracted beams interfere in the object region to generate a unique modulation intensity pattern 

for each point in the (x, z ) plane. 
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CHIRPT signals are formed by collecting all of the object fluorescent light on a single element detector that 
integrates the emitted object intensity. The total signal is a superposition of temporal intensity modulation patterns 
with an amplitude determined by the local fluorophore concentration, C(x,z), and is expressed mathematically as: 
S(t) = f dzdx Im(x, z;t) C(x, z) . Assuming the emitted fluorescence intensity has the same temporal structure 
as the illumination pattern, the instantaneous frequency of the signal recorded from the photodetector encodes the 
location of a fluorescent emitter in the (x, z) plane. The instantaneous frequency for a single emitter is: 


v(x, z ; t) 


d ( A$(x, z: t) 1 
dt \ 2ir J 

Ax M - 


d f AOPL(x, z; t ) \ 
= dt { A J 

ri2 A A z ( M Ki) 2 t 

71-1 \J 1 — (A M n\ t/ni ) 2 


1 dA$ a (t) 
2tt dt 


(3) 


The average instantaneous frequency of the illumination modulation encodes the transverse emitter position, Ax, 
while deviations in the instantaneous frequency over the scan time provide a unique determination of the axial position, 
Az. 

For an object composed of many emitters, the temporal patterns from each emitter in the object add together, 
and they must be separated to recover the object spatial information, C(x,z). Fluorescent emitters that lie in the 
object image plane, Az = 0, exhibit a uniform modulation frequency that depends on the lateral position, Ax. 
The temporal pattern for emitters axially displaced from the object plane carries an additional temporal phase. In 
CHIRPT, the temporal domain is proportional to lateral spatial frequency, so the temporal CHIRPT signal represents 
a superposition of frequency-domain Fresnel zone plates (FZP) in the paraxial optics limit. By noting that holography 
forms images through superposition of spatial FZPs im we can use an algorithm similar to that used for holographic 
imaging to numerically construct CHIRPT images (Supplementary Information). The CHIRPT algorithm selects the 
positive frequency sideband of the signal, S(t), then applies a propagation spatial frequency phase to recover the 
object, C(x, z). 

CHIRPT transfers the phase evolution of coherent illumination beams propagating through the object region onto 
a temporal pattern of fluorescent light emission. The spatial phase information can be recovered by numerically 
processing the CHIRPT signal, allowing Cl techniques to be directly applied to incoherent light. The CHIRPT 
signals can be numerically refocused to reconstruct object information in the entire (x, z) plane. The phase information 
also directly records imaging system aberration phase that can be numerically removed, eliminating optical imaging 
aberration distortions from CHIRPT images (Fig. S4). CHIRPT signals simultaneously collect fluorescent emission 
from a large range of depths in a single measurement, allowing the collection of out of focus contributions that would 
not be captured with a conventional imaging system, as well as digital refocusing of blurry objects. The result is that 
CHIRPT images have a large depth of field. 


RESULTS 


For a detailed description of the experimental methods, please consult the Supplementary Information. Briefly, we 
used a continuous-wave laser with a 532 nm wavelength to illuminate various fluorescent objects. Modulation masks 
with density AA: = 70/mm were used for all images shown in this work. The magnification and numerical aperture of 
the imaging system varied with the chosen objective lens. In the data presented here we used either a 50x/0.8 NA 
or a 20x/0.45 NA objective lens, resulting in an overall system magnification of 95x and 30.4x respectively. Both 
objectives were air-immersion, so in the data presented here we set n\ = = 1. Fluorescent light was measured in the 

epi-direction with a photomultiplier tube (PMT). The signal from the PMT was electronically amplified and filtered 
prior to digitization with a data acquisition (DAQ) card. All components of the microscope and data collection were 
controlled with a custom C=ff= application written in house. Please consult Table. SI for a summary of data acquisition 
times for each image presented here. 

To validate the CHIRPT theory, we imaged a single 100-nm-diameter fluorescent nanodiamond (FND) in the (x, z) 
plane. The fluorescence intensity distribution through the focus was measured by axially translating the FND along 
the direction of light propagation and recording the CHIRPT signal at each location. The complex CHIRPT signal 
was obtained from the measured voltage by a Hilbert transform (Fig. S3). The amplitude and phase of the complex 
CHIRPT signal are shown in Fig. s* and Fig. s> respectively. Equation [T] predicts that the phase evolution as a 
function of defocus should follow the Ewald circle shifted by the axial phase of the undiffracted beam. To validate this 
prediction, we examined the phase difference for two defocus positions relative to the nominal focal plane, indicated 
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by dashed lines in Fig. s* and Fig. n>. The phase differences for the positive and negative defocus values are shown 
in Fig. E with solid lines, while the theoretical phase difference predicted by Eq. [T] is plotted in dashed lines. The 
strong agreement between the measured and theoretical phase differences verified that the defocus phase was encoded 
in the temporal measurement as predicted by the theory. 
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FIG. 3: The spatial phase difference between the illumination beams is encoded in the CHIRPT signal, (a) 
Amplitude and (b) phase of the complex CHIRPT data from a fluorescent nanodiamond as a function of physical 
defocus, (c) Positive and negative defocus phase with respect to the phase measured at the nominal focus (blue and 
red solid lines respectively). Shaded areas denote the standard deviation obtained by measuring 20 consecutive 
scans. Dashed lines correspond to the theoretical Ewald phase for the defocus values read from the axial translation 
stage encoder. Data was collected at a scan rate of 17 frames/s. The total acquisition time for all 20 images was 

119 s. 


Measuring the intensity of a sub-diffraction-limited fluorescent object as a function of system defocus gives the 
point-spread function (PSF) and the corresponding optical transfer function (OTF) for our CHIRPT microscope in 
the [x, z) plane. Both the PSF and OTF are displayed in Fig. [4] The curvature of the OTF is further confirmation 
that the axial phase evolution in the CHIRPT microscope is well described by the Ewald phase. 



FIG. 4: CHIRPT encodes the Ewald phase as a function of physical defocus, (a) Fluorescence intensity vs. 
physical defocus measured by translating a 100 nm FND along the axial dimension. Data was collected with a 
95x/0.8 NA CHIRPT system. This image is the average of 25 images with a total collection time of 148.9 s. (b) A 
two-dimensional FFT of the image in (a) reveals the axial phase encoded as a function of defocus. The phase closely 
follows the Ewald phase (solid line). The Fresnel phase (dashed line) is also displayed to show that the paraxial 
approximation breaks down under these imaging conditions, so data must be numerically propagated with the 

Ewald phase. 
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To further confirm the spatial phase measurement from fluorescence intensity, we numerically reconstructed in-focus 
and out-of-focus 2D images in the (cc, y) plane of a single shell-stained fluorescent microsphere (Fig. [5]). The measured 
images can be propagated along the defocus distance by applying the conjugated Ewald phase, which is proportional 
to defocus, A z (Supplementary Information). The fluorescent sphere was first imaged in focus (Fig. [5^i), followed 
by an image collected after translating the bead in the defocus direction by -10.17 /mi, i.e., toward the objective 
lens. Images shown in Fig. B and Fig. B were computed by refocusing the collected images in Fig. B and Fig. B 
by -10.17 fim and +10.17 /tm respectively. Ideally, Fig B and Fig. B would be identical, as would Fig. B and 
Fig. [5}:. The excellent agreement between these images is further confirmation that the physical defocus phase is 
indeed encoded into the measured CHIRPT signal. Higher-order spatial phase distortions imparted by the imaging 
system were also recovered from other CHIRPT data and used to perform digital aberration correction (Fig. S4). 


Amplitude (arb.) 
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FIG. 5: Computed reconstruction of a fluorescent bead, (a) Fluorescence image of an in-focus 15 /tm shell-stained 
fluorescent bead, (b) Image of the same bead with -10.17 /tm of physical defocus (towards the objective lens), (c) 
Image formed by numerically propagating the image in (a) by -10.17 /tm in the axial dimension, and similarly, (d) 
the image formed by numerically propagating the image in (b) by +10.17 /tm. Images in (a) and (b) are averages of 
15 images, acquired at a frame rate of 17 frames/s with total acquisition time of 89.3 s. 


To explore the ability of CHIRPT to refocus fluorescent images in biological specimens, epi-fluorescent images from 
16-/tm-thick fixed murine intestinal tissue were collected. Figure [6|i shows a 2D image reconstruction formed by 
digitally refocusing a measured CHIRPT image by -4.5 /tm, while Fig. B shows the measured image. Both images 
show features, marked by arrows, being brought into focus at two different depths in the tissue by digital propagation 
of fluorescence intensity. To benchmark CHIRPT imaging in the tissue slices, the same region was imaged with a 
commercial spinning-disk confocal imaging system (Fig. [bj: and Fig. §1). The CHIRPT images clearly show features 
in focus at different axial planes, consistent with the confocal images. 

Comparison of the images in Fig. [6] reveals differences in the nature of image formation between CHIRPT and 
confocal microscopy. While a single CHIRPT image reveals cells from a range of depths in the tissue, a confocal image 
contains data from only a single axial plane due to optical sectioning. This numeric refocusing ability gives CHIRPT 
a very large depth of field (DOF). We computed the DOF in CHIRPT imaging using geometric considerations, and 
found it to vary inversely with the numeric aperture of the imaging system, i.e., DOF oc NA -1 (Supplementary 
Information). Thus CHIRPT is able to image the location of fluorophores in the object of a much larger axial range 
than is possible with conventional imaging, for which DOF oc NA -2 . 









FIG. 6: CHIRPT and confocal imaging of fluorescently-labelled murine intestine slices, (a) and (b) Digitally 
refocused CHIRPT images captured with 532 nm illumination light. Images were formed by propagating the average 
of 25 measured images. Panel (a) displays the image propagated by -4.5 ^m while panel (b) displays the measured 
image, (c) and (d) Confocal images with 561 nm illumination light. Images were collected with an axial separation 
of 4.5 /mi. In all four images, features that are in focus are denoted by arrows. Scale bars: 20 /mi. CHIRPT data 
was acquired at a rate of 17 frames/s for a total acquisition time of 443.7 s. 


We measured the DOF in CHIRPT by collecting an image of the illumination laser transmitted through a 20 
line/mm Ronchi ruling as the ruling was defocused (Fig.[7|. A 0.45 NA objective lens was used to resolve the object. 
The working distance of the objective lens is 510 pm, yet images were reliably formed though a scan range of 14 mm, 
or >27 times the working distance of the lens. Note that far away from the focal plane of the lens there is a loss of 
lateral spatial frequency content, made clear in Fig. [7]i by the odd harmonics of the binary, 50% duty cycle Ronchi 
ruling have a depth-dependent amplitude. The yellow line indicates the on-axis depth of field as a function of lateral 
spatial frequency, D0F 2 = 2w/(X f x ). The beam size, w, was estimated from the image in panel Fig. 0 to be 100 pm. 
Since the largest spatial frequency passed by the system is related to the numerical aperture by NA = A f x> ma x, the 
on-axis DOF for the full resolution of the lens can be recast as: DOF z = 2w/NA. 

We calculated the axial depth over which the full lateral spatial frequency support of the imaging system was 
retained to be approximately 440 pm, in good agreement with the data in Fig. [7] This is approximately 83 x larger 
than the predicted value of 5.3 pm for a conventional imaging system, defined as DOF ss 2nA/NA 2 |28j . 

For the majority of data presented here, we chose to display averaged images to improve the signal-to-noise ratio 
(SNR). This also helped to verify the theoretical analysis of the CHIRPT microscope, i.e., Fig. [3] and Fig. [4] We note 
however that CHIRPT reliably formed images with single-shot data acquisition. An example of this can be found in 
Fig. § While Fig. Hh shows a 2D image in the (x, z) plane reconstructed from data collected in a single rotation of the 
modulation mask, Fig. [8jo shows the average of 20 measured ID images. The single shot image was collected in 0.07 s 
while the averaged image was acquired in 141.4 s, not including time added by stop-and-hold stage scanning. While 
the SNR of Fig. [8 }d is indeed better than the single-shot image in Fig.[8ji, the single-shot image was acquired over 2,000 
times faster. The fidelity of the reconstructed image is remarkable given the relatively fast frame rate of the acquired 
data. The SNR is the primary factor limiting image acquisition speed in all forms of biological microscopy, and 
CHIRPT is no exception. However, CHIRPT shows great promise for drastically improving 2D frame rates. Another 
example of high-speed single-shot imaging with CHIRPT is shown in Fig. S6, where we imaged a fluorescent object 
with an acquisition time of 2.48 ms for a 2D (x, z) plane. The full 3D image volume was acquired in approximately 
one second, ignoring the hold time of the linear positioning stages. 

Finally, we note that the spatial resolution of CHIRPT imaging in the x dimension is diffraction limited. The 
maximal lateral spatial resolution of an imaging system, Sx, is limited by the largest lateral spatial frequency that 
can be passed by the lens, f x<max . The NA of the objective lens is directly proportional to the maximal lateral 
spatial frequency that the lens can accommodate. To achieve the maximal resolution afforded by the objective lens 
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FIG. 7: Large depth of field (DOF) imaging of a 20 line/mm Ronchi ruling with CHIRPT. (a) Transmissive image 
as a function of depth, where a depth of 0 mm nominally corresponds to physical contact of the ruling with the 
illumination objective lens, (b) The ID FFT of the image in (a) with respect to x, showing the decay of harmonics 
of the Ronchi ruling’s fundamental spatial frequency with defocus. Data was measured at 31 frames/s, for a total 
image acquisition time of 3.23 s for the single-shot image shown in (a). 


in an imaging configuration, f Ximax must contribute to image formation, which occurs when the back aperture of the 
objective lens is filled. Since CHIRPT images are formed by illumination with varying lateral spatial frequencies as 
a function of scan time, which are generated by the spinning modulation mask, optimal lateral spatial resolution is 
achieved when the highest spatial frequency generated by the modulation mask and image relay optics is greater than 
/a,max- In that case, Sx is determined by the NA of the imaging objective and the wavelength of the illumination 
light, A. Using the Abbe criterion, the resolution limit is Sx = A/(2NA). 

Conversely, we have constructed CHIRPT imaging systems for which the back-aperture of the objective lens is not 
filled. In that case the effective NA of the imaging system is limited by the largest spatial frequency generated by 
the modulation mask and the magnification of the imaging system. The effective NA of the system is related to the 
magnitude of the maximal lateral spatial frequency in the object region via: NA = A f x ,ma.x = ^2 sin # 2 (±T/2), and 
the lateral resolution of the system is: 


, Til 1 
Sx = ————r- 
n-2 M A k 


(4) 


For all data presented in the manuscript, the back aperture of the objective lens was filled - allowing imaging with 
diffraction-limited spatial resolution determined by the NA of the objective lens. The optimal lateral spatial resolution 
of each imaging system was computed using the Abbe criteria, results in a computed spatial resolution limit of 591 nm 
for the 20x/0.45 NA lens, and 333 nm for the 50x/0.8 NA lens. Note that while these resolution numbers correspond 
to uniform lateral spatial frequency support, the practical limit of the spatial resolution in CHIRPT is the amplitude 
of the measured signal as a function of scan time (spatial frequency). 

In the y direction, the spatial resolution is set by the effective NA determined by the filling of the pupil aperture 
along the y direction in the objective pupil plane. The y-resolution limit is also 5y = A/(2NA). 
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FIG. 8: Single-shot vs. averaged image acquisition of fluorescence form an FND. (a) A 2D image reconstructed 
from a single measured ID image at 14.3 Hz. (b) An image of the same region in the (x, z) plane formed by 
physically defocusing the FND and averaging twenty ID images at each defocus location. Total acquisition times for 

(a) and (b) were 0.07 s and 141.4 s respectively. 


DISCUSSION 


CHIRPT introduces a fundamentally new approach to imaging: the transfer of coherent propagation phase of 
illuminating beams to incoherent light emission. This phase transfer enables use of a powerful set of Cl tools that 
rely on spatial phase information to rapidly and efficiently record object information. These methods were previously 
not applicable to incoherent light imaging, which encompasses many of the most powerful tools in biological imaging, 
such as fluorescence and Raman scattering. We have demonstrated this new method and its ability to capture 
coherent phase information with recorded incoherent light using holographic imaging, which relies critically on phase 
information. 

The approach in CHIRPT can be applied to a broad range of imaging modalities that use waves to form im¬ 
ages. CHIRPT is highly achromatic, and can be applied to imaging radiation ranging from THz waves to x-rays. 
CHIRPT could prove particularly valuable for mid infrared (MIR) and terahertz (THz) imaging, where there is a 
dearth of imaging detectors. Moreover, the fundamental principles can be adapted beyond electromagnetic radiation, 
including acoustic waves used in seismic and ultrasound imaging, or in the hybrid acoustic/optical imaging technique 
photoacoustic microscopy m- 

CHIRPT has particular value for imaging in biological specimens. We have demonstrated high speed imaging of 
fluorescent and absorptive objects, with the ability to collect fluorescent images in a backscattered (epi) direction. 
Epi-fluorescent collection in CHIRPT has the potential to greatly impact biology studies, as it allows for light-sheet- 
microscopy-like imaging |29l I30| with a single lens. Furthermore, the image is acquired with a single-pixel detector, 
which allows for robust imaging in the presence of optical scattering, such as with biological tissues. The single 
pixel detector, in combination with high-speed modulation allows for transformative improvements in imaging speed. 
We have shown (x,z) plane collection in an epi-fluorescent configuration at 400 planes/second (Fig. S6), and have 
achieved >600 planes/second with multiplexed modulation masks (data not shown). CHIRPT allows collection of 
highly defocused light, enabling imaging volumes much larger than conventional approaches. An exceptional depth of 
field, approaching 30 times the working distance of the objective lens was demonstrated - well beyond that capable in 
conventional fluorescent imaging methods, particularly for 3D imaging. This set of capabilities brings new possibilities 
for multi-scale spatial and temporal imaging of biological organisms. 

Application of this technology to high speed imaging of biological dynamics and to large-scale imaging of biological 
specimens, including large tissue regions and ultimately organs, will provide a new window that could help form a 
comprehensive understanding of biological behaviors. Not only does this advancement enable the powerful toolbox of 
Cl techniques to be applied to imaging with incoherent light for the first time; it also enables the practical capabilities 
described above that permit new abilities in imaging. 
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SUPPLEMENTARY INFORMATION 

EXPERIMENTAL METHODS 


CHIRPT Microscope 


Two different continuous wave diode-pumped solid state (DPSS) lasers were used in these experiments, both of which 
operate at a nominal wavelength of 532 nm. Data in Fig. 7 and Fig. ^5]were acquired with a 100 mW source (Coherent 
Inc., Compass 315M-100). All other data was collected with a Sprout 5G pump laser (Lighthouse Photonics). 

Each laser was spatially filtered with a 10x/0.25 NA objective lens (Zeiss A-Plan, UIS) and a 10-pm pinhole 
(ThorLabs, P10S). The output of the spatial filter was collimated with an 80-mm focal length achromatic lens 
(ThorLabs, ACH254-080-A-ML), yielding a beam with a 4.95 mm full-width, measured as the width between the 20% 
and 80% peak intensity points. The remainder of the microscope was different for each laser source. 

For data collected with the Compass laser, the spatially-filtered beam was focused with a 300-mm focal length 
cylindrical lens, creating a line focus that extended laterally across the modulation mask (Fig. [l|, with a vertical 
height (full-width at the 20%-80% measurement points) of 15 pm. After passing through the modulation mask, the 
modulated line focused beam was re-imaged and demagnified with a 250-mm tube lens (ThorLabs ACH-254-250-A- 
ML), and an objective lens. The objective lens used for collecting experimental data depended on the particular 
experiment. Data in Fig. 7 was collected with a 20x/0.45 NA Zeiss A-Plan UIS objective lens, resulting in a system 
magnification of 30.4x. High speed data presented in Fig. was collected with a 10x/0.25 NA Zeiss A-Plan UIS 
objective lens, providing a system magnification of 15.2x. 

To increase the NA of the CHIRPT imaging system, the microscope was redesigned and rebuilt to allow for more 
precise control and alignment of the imaging system. In the redesigned CHIRPT system, we used the Lighthouse 
Photonics laser, as it has a relative intensity noise (RIN) significantly lower than the Compass laser. This substantially 
reduced the background noise on our CHIRPT data. This configuration was used to collect data displayed in all figures 
aside from those noted above. 

Since the exit pupil of high NA objective lenses often lies inside the objective case, access to the image plane 
conjugate to the object region (and modulator mask), where the horizontal slit was used to spatially filter the 
diffracted orders from the modulator (Fig. 2a), was not accessible in the initial CHIRPT microscope design. To 
overcome this limitation we designed the microscope with dual image relay systems to allow unrestricted access to the 
conjugate image plane. After spatial filtering and beam expansion, the beam was brought to a horizontal line focus on 
the modulation mask with a 150 mm focal length achromatic cylindrical lens (ThorLabs, ACY254-150-A-ML). The 
first image relay system consisted of a 125 mm focal length lens (ThorLabs, ACH-254-125-ML) after the modulator 
mask, followed downstream by a 100 mm focal length lens (ThorLabs, ACH-254-100-ML). The horizontal slit used to 
reject one of the diffracted orders was placed in the plane conjugate to the mask plane, which was located between 
these two lenses. The spatially filtered intensity profile was then image relayed to the object region with a 250 mm 
focal length tube lens (ThorLabs, ACH-254-250-ML) and an objective lens. For Fig. 3-Fig. 6, Fig. 8, Fig. and 

Fig. ^4] a 50x/0.8 NA objective lens was used (Zeiss, N-Achroplan 50x/0.8 NA Pol), and the overall magnification 
of the system from the modulator mask to the object plane was 95 x. 
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Modulation Mask 


Modulation masks were generated with a rotationally-varied groove density according to the relation: m(R, ip) = 
1/2 {1 + sgn[cos(Afc R <p)]}, where ( R , ip) represent polar coordinates relative to the center of the modulation pattern. 
A rendering of the mask is shown in Fig. The modulation pattern was printed in aluminum onto a glass substrate 
(Projection Technologies, Dallas, Texas) at 3600 DPI. A center hole was cut with a water jet so the disk could be 
mounted on a custom chuck, and attached to a motor (Faulhaber, 2057S012BK1155). The desired angular velocity of 
the motor was set by the user in the data acquisition software and controlled by an external speed controller module 
(Faulhaber, MCBL3006). 


774 



FIG. SI: A CHIRPT modulation mask with Afc = 2/mm. As the mask is spun at constant angular velocity 
uj r = 2nu r , the modulation frequency varies linearly as a function of radial coordinate, R. Since the line focus is 
aligned along the x axis, the modulation frequency experienced by the illumination light sheet is modulated with a 
temporal frequency that varies linearly with lateral position. Temporal modulation frequency, i/ t , and lateral 
position, x, are related by the density of the mask, Afc, and the rotational frequency, v r , by the expression: 
v t = k-[ x = Afc v r x, where K\ = Afc v r is known as the chirp parameter. The inner and outer radii of the 
modulation pattern are 9.5 mm and 31 mm respectively. The green rectangle on the right side of the mask indicates 
the location of the laser line focus incident on the modulation pattern. 


To measure and correct for small modulation frequency errors arising from modulation mask centering misalignment, 
which inevitably occurs during mounting of the modulation disk, the illumination laser was brought to a focal point 
on the modulation mask with a second cylindrical lens with a focal length of 125 mm (Thorlabs, ACY254-125-A-ML) 
that was oriented perpendicularly to the permanent cylindrical lens. The throughput laser intensity as a function of 
scan time was measured with a photodiode for 512 traces. For each trace, a spectrogram was computed and the central 
modulation frequency was obtained by fitting the centroid of the spectrogram (Fig. ^2j. The average modulation 
frequency was then computed from the set of 512 centroid frequency traces. If the modulation mask was mounted 
perfectly, the modulation frequency would be constant as a function of scan time since the focal spot samples the 
same position in space. Deviations from constant frequency correspond to imperfect mounting of the mask on the 
rotation axis of the motor. Mathematically, we can treat these deviations as a time dependent phase shift, such that 
the signal measured with the photodiode has the form: 

V m (t) oc cos [27 t t x s + <j) m (t)\ = cos [2 tt y s t + 4>m{t)\ (5) 

where K\ is the chirp parameter of the mask, x s is the lateral location of the focal spot on the mask, and <j> m (t) is 
the phase that causes deviations from a constant modulation frequency, v s . This phase can be recovered from the 
modulation frequency measurement by numerically integrating the centroid frequency vs. time data and multiplying 
by 27 t (cf. Eq. (2) of the manuscript). The conjugate of the mask phase was applied to measured CHIRPT data to 
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correct for the phase disturbance. 



Scan time (ms) x (pm) 

FIG. S2: Disk phase aberration correction, (a) A spectrogram of the modulation frequency of the laser intensity 
at a single lateral location as a function of scan time. The computed centroid of the modulation frequency is shown 
by the overlaid red line, (b) The mean modulation frequency (solid line) and standard deviation (shaded region) 
were computed from 512 temporal scans, (c) The phase imparted by the modulation mask was uncovered by 
numeric integration of the mean modulation frequency, (d) A measured {x, z ) image of a 100 nm FND without 
correcting for the disk aberration phase, (e) A corrected image, computed by applying the conjugate disk aberration 
phase to the data, (f) The intensity of the corrected (solid line) and uncorrected (dashed line) fluorescence intensity 

near the focal plane. 


All data presented in this work used mask designs with Afc = 70/mm. For high-speed data in 
four masks was utilized to increase imaging speeds. 


Fig. f] a disk with 


Detection System 

Fluorescent light intensity was measured with a photomultiplier tube (PMT). Data in Fig. f] were recorded with 
a Hamamatsu H9305-03. All other fluorescent images were measured with a Hamamatsu H7422P-40. 

For the high-speed data in Fig. ^5] the signal from the PMT was high-pass filtered by a current preamplifier 
(Stanford Research Systems, SR570) with 10 kHz corner frequency and 12 dB/octave rolloff. Data in Fig. 3-Fig. 6, 
Fig. 8, Fig. f] and Fig. €1 were collected with a low-noise current amplifier (DHPCA-100, Femto Messtechnik 
GmbH, Berlin, Germany). Transmitted laser intensity presented in Fig. 7 was measured with a large-area photodiode 
(ThorLabs, DET100A). Amplified and filtered signals were digitized with a PCI data acquisition card (National 
Instruments, PCI-6110 with BNC-2090A breakout box). Data was displayed to the user and saved for post-processing 
with a custom C# Windows application. 


Fluorescent Objects 


A 100 nm diameter fluorescent nanodiamond (FND) was used as a fluorescent probe for data presented in Fig. 3, 
Fig. 4, Fig. 8, Fig. ^3] and Fig. £|4j A 10 /xL drop of FND slurry (ND-400NV-100nm-10mL, Adamas Nanotechnologies, 
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FIG. S3: CHIRPT image reconstruction, (a) A single temporal scan collected by measuring fluorescent light from 
a fluorescent nanodiamond (FND) near the focal plane of the microscope with a 50x/0.8 NA objective lens, (b) The 
spectral density of the measured temporal trace, (c) An inverse FFT (IFFT) is used to recover the complex 
CHIRPT trace in the temporal domain, Sj + (t). (d) The CHIRPT trace is demodulated by the carrier frequency to 
remove the linear phase. Other phase disturbances, i.e., systematic aberration phase, or pupil phase, and the disk 
aberration phase are also removed to compensate for known aberrations to the system, (e) The demodulated and 
phase-corrected CHIRPT data is digitally propagated by applying the Ewald phase, (f) A one-dimensional FFT 
with respect to the lateral spatial frequency domain reveals a numerically-reconstructed two-dimensional image in 
(x, z ) space, (g) An image in the (x, z) plane formed by physically scanning the FND along the defocus dimension. 


Raleigh, NC) was drop cast onto a standard microscope slide and allowed to dry at room temperature. A cover slip 
was then placed over the sample to minimize spherical aberration. 

A single 15-/rm-diameter, shell-stained fluorescent microsphere in a prepared slide (LifeTechnologies, FocalCheck™ 
Slide 1, Well Al) was used for CHIRPT data presented in Fig. 5. Although the slide contains fluorescent microspheres 
with a variety of fluorophores, the orange fluorescent spheres were most visible with the illumination wavelength of 
532 nm. 

The murine intestinal images presented in Fig. 6 were collected from a commercial prepared slide (LifeTechnologies, 
FluoCells®Prepared Slide =^4). 


CHIRPT Reconstruction 


Data processing was performed in MATLAB (The MathWorks, Natick, NJ). Image recovery for a single temporal 
data trace was achieved by a fast Fourier transform (FFT) operation. The lateral spatial axis was calibrated according 
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to the relation f x x = v t t, which yields: 


n i v t n i v t n 1 v t 

x = — —— = — -=- (o) 

ri2 M Kl 712 M A K I/ r U2 K2 

where is the temporal modulation frequency, v r is the rotation frequency of the mask, A k is a known design 
parameter of the mask, and M is the magnification of the illumination microscope. The chirp parameter in the object 
region, k 2 = M was determined empirically by collecting CHIRPT signals while translating an isolated fluorescent 
object laterally, i.e., in the x direction. The fluorescent object was typically a 15 /iin shell-stained bead. The centroid 
modulation frequency of the fluorescence distribution was computed for each CHIRPT trace, which corresponded to 
a known lateral location. The chirp parameter was determined from the slope of a linear fit to the centroid frequency 
vs. lateral position data. 

Refocusing of fluorescent light intensity was performed by applying a propagator in the lateral spatial frequency 
domain, which is directly proportional to time. The reconstruction processes is summarized visually in Fig. £|3] 
Propagation was applied to the positive frequency complex sideband of the measured CHIRPT time trace. This 
sideband was selected by filtering out only the positive temporal Fourier sideband that contains the modulation 
frequencies. The centroid of the positive frequency sideband was computed and used to demodulate the temporally- 
modulated complex data. The demodulated data were numerically propagated with a propagation phase conjugate 
to the axial phase component in Eq. (1) This phase varies linearly with defocus distance, A z. After application of 
the Ewald propagation phase, a ID FFT with respect to lateral spatial frequency then gives the propagated image. 
Mathematically, propagation can be described as: 


S[+(x,z) =J r f x <( 5i + (/ a ,, 2 )exp 


=Fi 27 t ^ Az 
A 




(7) 


where the superscript (r) is used to denote the refocused data, is the Fourier transform operator with respect to 
lateral spatial frequency, and Si+(f x ) denotes the complex CHIRPT data, which is synthesized from the measured 
real-valued data in the temporal domain, S±(t), by a Hilbert transform. The sign of the argument used to propagate 
depended on which diffracted order was passed through the spatial filter (i.e., j = ±1). 

Data presented in Fig. ^5] was further processed by multitaper power spectral density (PSD) analysis using the 
function pmtm() in MATLAB with a time-bandwidth product of 5/2. The processed image was computed by taking 
the square-root of the PSD estimate. 


TABLE I: Data Acquisition Parameters 


Figure 

v r (s 1 

) N scans 

Ap ix 

^total (s) 

V 

3 

17 

20 

101 

119 

2D 

4 

17 

25 

101 

148.9 

2D 

5 

17 

15 

101 

89.3 

3D 

6 

17 

15 

101 

443.7 

3D 

7 

31 

1 

100 

3.23 

2D 

8a 

14.3 

1 

1 

0.07 

2D 

8b 

14.3 

20 

101 

141.4 

2D 

S4 

14.3 

20 

101 

141.4 

2D 

S5 

403 

1 

401 

0.995 

3D 


Most data displayed in the main text were formed by averaging multiple images. Table ^summarizes the acquisition 
parameters for each image in this work, where v r is the frame rate of ID data acquisition (2D after CHIRPT 
reconstruction), N scans is the number of ID data scans averaged for each position that the object was translated, 
Api x is the number of positions at which data was collected, i t ot is the total image acquisition time, and T> is the 
dimensionality of the reconstructed CHIRPT data. Note that the same data is displayed in Fig. 8 and Fig. ^3] 
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Confocal Imaging 


Spinning-disk confocal images of murine intestinal tissue were collected with a commercial microscope base (1X81, 
Olympus) equipped with a spinning disk confocal head (CSU22, Yokogawa Instruments). Fluorescence was excited 
with an illuminating laser with 561 nm wavelength. An objective lens with NA of 0.5 focused the illuminating light 
into the specimen (Olympus UPlanFL N 20x/0.5 Phi). Fluorescent light was isolated from the illumination laser with 
an interference filter that had a 73 nm bandwidth centered at 617 nm (FF02-617/73, Semrock). Fluorescent images 
were formed on an electron multiplying charge coupled device (EMCCD) camera (Cascade 11:1024, Photometries). 
Image acquisition was controlled with SlideBook v.6 software (Intelligent Imaging Innovations). 


THEORETICAL ANALYSIS OF CHIRPT 


In CHIRPT, the spatial phase of coherent illumination beams is imprinted onto a temporal intensity modulation 
pattern that is unique to each (a:, z) position. Here we compute the illumination intensity pattern by assuming 
spatially-coherent, monochromatic plane waves interfering in the object region. We then show that the spatial phase 
difference between the two interfering illumination beams is encoded as temporal modulations on the illumination 
intensity pattern. 

CHIRPT encodes complex spatial image information into a temporal intensity pattern by measuring the radiant flux 
from the contrast distribution on a single-element photodetector as the illumination intensity pattern is modulated 
over time. The CHIRPT signal has the form 

/ OO 

d 3 r Ain(r; t) C(r, t) (8) 

-OO 

where /jh (r: t) is the illumination pattern in the object region, and C(r,t) is the contrast distribution. The contrast 
distribution describes the method by which illumination intensity is transferred to contrast intensity. For example, 
the contrast distribution can simply represent the transmittance of the object if the measured contrast intensity is 
the illumination source after propagation through an object displaying absorption. In that case, C(r,t) = T(r,f) = 
1 — A(r, t), where T is the intensity transmission of the object, and A is the intensity absorption of the object, and an 
image is formed from the object transmission. The contrast can also represent the spatial concentration distribution 
of fluorophores in an object, scaled by the absorption cross section for the excitation light and efficiency of fluorescent 
emission. Throughout this analysis, we shall assume that the contrast function is constant for the duration of a 
scan, so that C(r,t) —> C{ r). Moreover, we will assume that the contrast function is linearly proportional to the 
illumination intensity. 

To compute the form of the illumination intensity, which arises from the spatial interference of the undiffracted 
and diffracted beams, consider two plane waves interfering in the lateral-axial plane, (x,z), as shown in the object 
region in Fig. la. We define the coordinate system such that the undiffracted beam (green) lies along the optic axis, 
z. Throughout the following analysis, we use a subscript pair (j, s) to denote the diffracted order j in region s, where 
s = 1 corresponds to the mask region and s = 2 corresponds to the object region. For example, the electric field of 
the undiffracted beam (j = 0) in the object region is: 


v 0 , 2 (x,z;t) = ao e lk °’ 2 ' r = a 0 e lk2Z (9) 

where a o is the amplitude of the field, ko .2 is the wave-vector of the undiffracted beam in the object region, = 
|k 2 | = 27rn2/A is the wavenumber, and n 2 is the refractive index in the object region. 

Conversely, the first-order diffracted beam (j = +1) propagates at an angle with respect to the optic axis that 
varies with scan time, 0 2 (i), so the wave-vector for the positive diffracted beam is: 

ki, 2 (t) = fc 2 sin e 2 (t)± + k 2 yl- sin 2 0 2 (t) z. (10) 

Additionally, the diffracted beam accumulates an optical path length (spatial phase) that varies with propagation 
angle, and corresponds to the pupil phase of the imaging system, & pU pii(fx)- Since the lateral spatial frequency, f x , 
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is directly proportional to scan time, t, the electric field of the diffracted beam in the object region is: 


Vi : 2 (x,z;t) = ai e lkl ’ 2 ' r 


= ai exp [i $ pU pii(t)] exp [i k 2 x sin 9 2 (t)] exp 



( 11 ) 


From Fig. la, it is apparent that the propagation angle in the object region, 9 2 (t), is related to the angle of diffraction 
from the mask, 9i(t). Specifically, these angles are related by the magnification of the image relay system, M, which 
is the ratio of the tube lens focal length to the objective lens focal length, i.e., M = F t /F 0 . Moreover, the diffraction 
angle in the mask region is related to the spatial frequency of the modulation pattern on the mask according to the 
expression: ni sin(?i(f) = \f x ,i(t), where is the spatial frequency of the mask, and ni is the refractive index 

of the medium that the mask is immersed in. Altogether: 

sin 9 2 {t) = sin 9 X (t) = — Mf x> i(t). (12) 

F 0 ni 

To derive the time-dependent spatial frequency, we consider the modulation mask. In polar coordinates, the 
modulation mask has the form: 

m{R,ip) = ^ + isgn[cos(Afci?^)] (13) 

where A k describes the density of features on the mask and has units of spatial frequency, and the function sgn[-] is 
included to account for the binary amplitude modulation provided by the modulation mask. This expression represents 
the full two-dimensional pattern of the mask, while we are interested in the region sampled by the line focus, located 
at y = 0, and approximately centered on one side of the mask such that the illumination distribution in x lies between 
the inner and outer radii of the mask (Fig. To account for the changing modulation pattern sampled by the line 
focus during mask rotation, and thus find the local spatial frequencies as a function of scan time, we consider the 
mask: 


m(R, <p - ip 0 ) 


\ + 7 jSgn{cos[Afc R (v? - y> 0 )]} 


(14) 


where (po describes the rotation angle of the mask. Since the mask rotates at a constant angular velocity oj r = 2 ti v r , 
where v r is the rotational frequency of the mask, the rotation angle is: ipo = 2 tt v r t. 

To uniquely determine the spatial phase of the illumination microscope, only the first diffracted order from the 
mask is permitted to interfere with the undiffracted beam in the object plane. Higher diffracted orders resulting from 
the binary modulation scheme are omitted in the following analysis, and also excluded in our CHIRPT experimental 
design, so that the mask we consider is: 


m(R,<p- (p 0 ) 


\ + \ cos[Afc R{<p- ¥> 0 )] = 1 + 1 C0S I V’o) ] 


(15) 


The linear spatial frequency on the mask in the lateral dimension is ED: 


r u\ 1 de />( T ^( t )) 
2n dx 


(16) 


The local region sampled by the beam is along the line where ip = 0, so the local linear spatial frequency imparted to 
the beam is: 




2n 


COS(f 


dr 


1 . dcj>{ v,<p) 

- smip — - - 

r o<p 


= Kit 


(17) 


- (fi =0 


where K\ = Afc v r is the so-called chirp parameter, which describes the linear increase in modulation frequency with 
lateral position on the mask. Equation (121 becomes: 

sin 0 2 (t) = — M A k v r t = — (M k± t). 

ni n\ 


( 18 ) 
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The lateral spatial frequency in the object region is then: 


fx, 2{t) 


ri 2 sin 6*2 (t) 
A 


— M Kit 
n i 


( 19 ) 


The electric field of the diffracted beam can now be written: 


Vi i 2 (x, z;f) = ai exp 


i 27 t ( — M K\t \ x 


exp 


i 27 r x*V 1- 


ni 


M K\t 


pi ^pupil (*) 


= ai exp[i^ latera i(a;;t)] exp[i </ axia i(,z; t)] exp[i $ pup ii(t)] 


( 20 ) 


Here we see that the diffracted order accumulates a linear spatial phase that varies with lateral position, a circular 
spatial phase that varies with axial position, and a pupil phase that varies only with scan time. 

The axial phase component describes the optical path length (OPL) change imparted to the diffracted beam as it 
is scanned across the pupil as a function of time. The magnitude of the linear spatial frequency in the object region 
is / 2 = k 2 /{ 27r) = n 2 /\, and the axial phase can be written as: 


^axial(^j 6) — k 2 Z 



(—MKit\ =27 rf 2 z\ 


'fxAty 

h 

w J V 



( 21 ) 


This phase corresponds to physical defocus in the imaging system, which is represented schematically by the curved 
orange wavefronts in Fig. la, and is consistent with the angular spectrum representation [32] • This phase is formally 
equivalent to the Ewalcl phase. 

We note that it is possible to construct an imaging system for which the paraxial approximation applies. Applying 
a small-angle approximation to this expression, wherein we assume sin 2 d 2 (t) = [(AM k\ t ) /ni ] 2 <C 1 , or equivalently, 
[fx,2 if) / /g] •C 1 , the axial phase can be written as a quadratic phase variation: 


axial (*£; A 6) 




_A 

n 2 


2 


fx, 2(6) 


( 22 ) 


This phase variation is equivalent to the Fresnel propagation phase eh, and corresponds to physical defocus in the 
imaging system, i.e., axial displacement away from the object focal plane. 

Finally, the illumination intensity in the object region is computed by a coherent summation of the diffracted and 
undiffracted beams: Ia\(x,z;t) = |uo j 2 (a:, z\ t) + Vi ) 2 (x, z; 6 )| 2 . We assume beam amplitudes ao = 1/2 and ai = 1 / 7 T 
based on the diffraction efficiency from a square grating. The illumination intensity evaluates to: 


Im(x,z;t) = I 0 + h(x,z]t) = 
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1 

+ — cos 
7 r 


2n 
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2n 
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n 1 


A<frpupii(t) 


(23) 


The illumination intensity consists of a constant (DC) term, Iq, and a temporally varying (AC) term, Ii(x,z\t). The 
CHIRPT image is encoded by the AC component of the illumination intensity: 


Si{t) 


da; d z Ii{x, z; t) C{x, z) 


(24) 


The argument of the cosine term in Eq. (231 contains a two-dimensional map of the spatial phase variation in the 
object region due to system defocus as well as the pupil phase of the imaging system. 


Digital Aberration Correction 


The differential phase between the two illumination beams encodes the pupil phase of the imaging system as well 
as physical defocus. The pupil phase accounts for deviations in the optical path length due to systematic aberrations 
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in the imaging system. Since this phase difference is encoded in the temporal modulation frequency of the measured 
signal, it can by digitally removed to correct these aberrations. 
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FIG. S4: Systematic aberration phase correction, (a) Measured phase in the lateral spatial frequency domain, 
collected with a 95x/0.8 NA CHIRPT microscope by physically defocusing a 100 nm FND. This phase map is the 
average of 20 scans, (b) Removing the theoretical Ewald phase from the 2D phase map in (a) reveals the residual 
phase as a function of defocus, (c) The average (solid line) and standard deviation (shaded region) of the residual 
phase in (b) representing the systematic aberration phase of the imaging system, (d) A corrected phase map, formed 
by removing the mean aberration phase in (c) from the measured phase in (a), (e) An image of the FND before 
removing the systematic aberration phase, and, (f) after removing the aberration phase, (g) The intensity of the 
uncorrected image (dashed line) and corrected image (solid line) at the location indicated by the solid black lines in 

(e) and (f). Total collection time for all 20 scans was 141.4 s. 


An example of systematic aberrations in an imaging system is shown in Fig. €1 where the fluorescence intensity 
emitted by a 100 nm fluorescent nanodiamond (FND) was measured as a function of system defocus. The average of 
20 phase measurements as a function of physical defocus was computed from the data. In the absence of pupil phase, 
the phase recovered from this measurement would be identical to the phase defined by the Ewald circle, i.e.: 


^Ewald (/)e,2 ; Z) = 2 tT — Az 



(25) 


To uncover the pupil phase we computed the Ewald phase, with Az = 0 at the nominal focus of the measured image, 
and removed it from the measured phase. The resulting phase map showed a relatively constant phase over a range of 
defocus values. The mean pupil phase was computed and removed from the measured data. Two-dimensional images 
of the FND fluorescence intensity in ( x , z)-space were computed from the measured data and the phase-corrected 


































20 


data. The symmetry of the intensity distribution vs. defocus was improved in the image when the residual pupil 
phase was removed. Intensities of each image near the focal plane show a significant improvement of the focal quality 
in the phase-corrected image. 


High-Speed CHIRPT with Fluorescent Molecules 

Since excited fluorescent molecules emit light spontaneously, there is a delay between excitation of a fluorophore 
and emission of a fluorescent photon. This delay is characterized by the fluorescent lifetime, Tf, of the molecule. A 
requirement of CHIRPT imaging is that emitted fluorescent light must follow the temporally-modulated illumination 
intensity pattern. This requirement is met when the modulation frequency of excitation is well below the inverse 
lifetime of the molecule, i.e., vt tJ 1 . If the temporal frequency of excitation does not fulfill this requirement, the 
emitted fluorescence intensity acquires an appreciable temporal phase shift, <fif, as well as a change in the modulation 
depth, rrif. The phase shift contributes to a phase delay in the measured signal, while loss of modulation depth 
reduces the signal-to-noise. 

The emitted fluorescence intensity is directly proportional to the excitation intensity, I ex (x, z;t), which is propor¬ 
tional to the product of the illumination intensity and the contrast distribution. Taking fluorescent lifetime into 
account, the temporally-varying fluorescence intensity emitted by a fluorescent molecule becomes: 


x , z\ t) oc m/ cos < 4>f + 2 tt Ax M t^j + k 112 z 


(1 _ ,'AMki* , _ 1 


n 1 


(26) 


where we have retained only the the portion of the fluorescence intensity that encodes the CHIRPT image - the 
temporally varying signal. 

Fluorescent lifetime imaging with temporally-modulated excitation sources has been studied elsewhere [55] . Fol¬ 
lowing this analysis, and consistent with the theoretical analysis of CHIRPT imaging presented above, we consider 
fluorophores excited by a sinusoidally-varying intensity pattern. Assuming single-exponential lifetime decays, the 
phase delay of the emitted fluorescence intensity in CHIRPT imaging evaluates to: 


and the depth of modulation is: 


4>f = 2lT KlX Tf , 


ITlf = \J 1 + (27T Kl X TfY 


(27) 


(28) 


Thus we can write the temporally-varying signal from the photodetector [Eq. (24)] as: 

OO 

S'i(t) oc JJdxdz \J1 + (27t M Hi x Tf) 2 

— OO 

x cos j 27r Kq x t/ + 27t Ax ^— M «q t^j + kn ,2 z ^1 


XM Kit 


- 1 


(29) 


The expression above indicates that for the recorded signal to reproduce the modulation frequency patten of the 
illumination intensity, the phase shift due to the fluorescent lifetime must be much less than unity: 2ttki xTf <C 1. 
Assuming a typical fluorescent lifetime of 10 ns [3D, we computed this phase delay for data presented in Fig. 
which had the highest temporal modulation frequency in this work, to be ~ 0.036 radians. 


Vertical Diffraction from the CHIRPT Mask 

The CHIRPT illumination intensity pattern is formed by allowing only one of the diffracted beams to interfere 
with the undiffracted beam in the object region. This is accomplished with a horizontally-oriented slit placed near 
the pupil plane of the illumination lens that passes only a portion of the undiffracted beam and a portion of one of 
the positive or negative diffracted orders (Fig. 2a). This separation is possible due to a non-zero spatial frequency 
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FIG. S5: High-speed epi-fluorescent CHIRPT images of lens tissue stained with highlighter ink. Scan times were 
measured by (a) using repetitive features in the temporal scan, and (b) adjusting the trigger delay of the DAQ card 
to collect two successive central peaks located near the center of the scan time, i.e., t = 0. Scan periods were found 
to be 2.48 ms in each case, corresponding to data acquisition for a full 2D (x, z ) plane at 403 frames/s. Each image 
is composed of 401 vertical scans for a total acquisition time of 995 ms for the full volume. In both (a) and (b) the 
temporal waveforms are the mean of 401 scans at varied vertical positions, (c) and (d), Digitally propagated images 
of the lens tissue, with the amplitude of each image estimated with a multi taper analysis (time-bandwidth product 
= 5/2). The propagation distances were -42 /jin and +20 /jm in panels (c) and (d) respectively. Scale bar: 100 pm. 


component in the vertical ( y ) direction that is independent of the rotation angle of the mask. Evaluating the local 
linear spatial frequency eh imparted by the modulator in the vertical dimension, we find: 


f _ 1 d</>{r,<p) 
fv>1 ~ 2n 8y 
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dr 
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(30) 


The non-zero spatial frequency in the vertical dimension causes the diffracted order to propagate at an angle with 
respect to the optic axis, i /q. In the pupil plane of the imaging system, the displacement of the diffracted beam from 
the optic axis is given by: 


F t sin^r = F t = F f 


n i 


A A k 

27T7ll 


(31) 


where F t is the focal length of the tube lens. This displacement permits insertion of the horizontal slit to select only 
one diffracted order. 


For the data shown in Fig. ^5] a modulation disk with 4 masks was used to increase image acquisition speed by 
providing a full temporal scan in only 90 degrees of mask rotation. To combine multiple mask patterns onto a single 
disk, the pattern was phase-wrapped in angle. Consequently, the density of features in the vertical dimension sampled 
by the incident line-focused beam increases as well. In general, the spatial frequency in the vertical dimension increases 
by a factor of N in this scenario, where N is an integer number of masks per disk. Therefore the displacement from 
the optic axis in the pupil plane is increased by a factor of A as well. 
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Depth of Field in CHIRPT 

One important advantage of CHIRPT imaging is an increased depth of field (DOF) as compared to conventional 
linear imaging. In tightly-focused imaging techniques, such as laser-scanning microscopy, the depth of field and 
lateral spatial resolution are both set by the numerical aperture (NA) of the illumination objective, and the DOF 
varies inversely as the square of the NA. CHIRPT does not suffer this same limitation because collimated beams with 
varied propagation angles are utilized to sequentially measure the spatial frequency content of an object. Since the 
CHIRPT signal is generated by the temporally varying part of the illumination intensity, Ii(x, z; t), the depth of focus 
is dictated by the region over which the intensity modulation is non-zero. This region is formed by the interference 
of the re-imaged diffracted orders. 

In the theoretical analysis thus far, we have assumed plane-wave illumination. This causes both the field of view 
(FOV) and depth of field (DOF) to be infinite. Let us assume instead that the illumination beam incident on the 
modulation mask has a Gaussian intensity profile in the lateral dimension and is collimated in the axial dimension. 
We may write the incident illumination field propagating along the optic axis as 

E- mc {x,z) = u{x)e ik ° z = e - {x - x °) 2/w2 e ikoZ (32) 

where w is the width of the Gaussian distribution. 

Due to diffraction from the modulation mask, the amplitude distribution propagates at an angle with respect to 
the optic axis in the object region (Fig. la). We denote the amplitude of a diffracted order j in the object region as 
Uj/(x , z; t). The amplitude of the undiffracted beam in the object region is simply given by the illumination amplitude 
scaled by the magnification of the illumination system: 

uq{x, z\t) = uq(x/M) = e -M 2 ^-x 0 ) 2 /w\ ( 33 ) 

The diffracted order, j = 1, propagates at an angle with respect to the optic axis that depends on the rotational 
position of the mask and therefore depends on time [Eq. ©]• To simplify the following computations, we apply 
the approximation that upon diffraction from the mask the amplitude distribution simply rotates by the diffraction 
angle, i.e., there are no distortions to the amplitude distribution. In the rotated coordinate basis (x',z'), where z' is 
the optic axis of the diffracted beam in the object region, the spatial distribution of the diffracted order is 

Ul (x'/M) = e ~ M2 (*'-*o) 2 /™ 2 (34) 

The DOF is dictated by the region of overlap of the diffracted and undiffracted spatial distributions: 

A(x, z\ t ) = uq(x, z; t) u±(x, z ; t) = uq(x) ui(x') (35) 

The diffracted order distribution can be written in the unrotated coordinate basis, (x,z), via the coordinate transfor¬ 
mation: 


x' = x cos 02 (t) — z sin 0 2 (t) 

Assigning the centroid of each distribution to the origin, xq = x' 0 = 0, we can write: 

x 2 cos 2 02 (t) z 2 sin 2 0 2 (t) 2 xz sin 0 2 (f) cos 02 (t) 
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Using the shorthand notation T = A M k\ t, the overlap distribution in the object region is: 


A(x, z, t) = exp 
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(37) 


(38) 


Although the DOF varies slightly with position in the object region, it is instructive to examine the DOF along the 
optic axis, i.e., x = 0: 


A (a; = 0, z, t ) = exp 


/ A M K\ t\ 2 2 


A » ) z \ 

= exp 


A f x ,2(t ) 


(39) 























23 


where f x , 2 (t) is the lateral spatial frequency in the object region. We define the DOF along the optic axis, DOF z , as 
the full-width of the overlap distribution at the 1/e-point, which we can extract from Eq. (391: 


DOF.(f) = 


2 w 


2 w 


AM Kit A / x>2 (t) 


(40) 


Since the highest spatial frequency that contributes to the measurement corresponds to the maximal scan times, 


t = ±T/ 2, Eq. (401 can be recast in terms of the numerical aperture of the imaging system to find the DOF for the 
full numerical aperture of the imaging system. 


DOF, = 


2 w 


4 w 


XMk i T/2 A M A k v r T 


4 w 

A M A k 


2 w 
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(41) 


Unlike conventional light focusing, where the DOF varies as 1/NA 2 , the DOF in CHIRPT at the full numerical 
aperture of the system varies as 1/NA. This allows for high resolution imaging from a much broader depth. 
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